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ORBIT DETERMINATION FOR LOW-THRUST SPACECRAFT; 
CONCEPTS AND ANALYSIS* 


Abstract 

This paper re-evaluates Earth-based orbit 
determination capability for SEP spacecraft in the 
light of recent developments both in multi •station 
tracking concepts and in thrust subsystem error 
modeling. Five different tracking strategies are 
applied to a 15 day segment of an Encke rendezvous 
mission. Both optimal and suboptimal orbit deter- 
mination performance are determined for a wide 
range of process noise .parameter values. The 
multi- station tracking techniques are found to be 
extremely effective, reducing orbit determination 
errors by orders of magnitude over that obtained 
with conventional single -station tracking. Explicitly 
differenced multi-station data (QVLBI) is found to 
be least sensitive to gross modeling errors, but if 
a reasonably good process noise model is available, 
explicit differencing is not required. 

L Introduction 

Solar electric propulsion (SEP) is characterized 
by high level stochastic nongravitational accelera- 
tions resulting from random variations in the thrust 
process. The random accelerations are roughly 
three orders of magnitude larger than those due to 
gas leaks, solar pressure variations, etc. on bal- 
listic missions and at such levels constitute the 
dominant error source for Earth-based interplane- 
tary navigation of SEP spacecraft. (1) 

Since reduction of thrust subsystem errors to 
ballistic levels does not appear to be feasible, the 
successful navigation of SEP missions depends upon 
making the orbit determination process more toler- 
ant of stochastic forces. As a first step toward 
this end it is clear that some form of stochastic 
model is required to represent the random thrust 
process, and this implies the need for a sequential 
filter (2) to accommodate the resulting process 
noise model. Unfortunately, even with optimal 
sequential filtering under the optimistic assumption 
that the process noise model is a perfect represen- 
tation of the physical random process, orbit deter- 
mination capability for SEP spacecraft remains 
grossly inferior to that for ballistic craft when con- 
ventional Earth-based tracking is assumed. 

Adaptive filtering techniques (^»4) may be useful 
for SEP missions as a means of identifying and 
updating the process noise model in real time. But 
it is clear that such techniques can only enable the 
filter to approach more closely in operation the 
theoretical optimal performance indicated when the 
process noise model is assumed correct and there- 
fore needs no updating. If the idealized (optimal) 
filter performance is unsatisfactory, as is gen- 
erally the case for SEP, no filter algorithm. 


however elegant, can remedy the situation. 
Something more fundamental is needed; i. e. , the 
information content of the data available to the fil- 
ter must be upgraded. 

One superficially attractive approach would be 
to reduce the uncertainties at their source, i. e. , 
to measure the actual nongravitational accelera- 
tions in real time with precision onboard acceler- 
ometers and make this data available to the orbit 
determination filter. Unfortunately, the instru- 
mentation requirements (resolution and alignment) 
for this approach to be effective are prohibitive. 
Specifically, to approach ballistic orbit determina- 
tion accuracy the acceleration uncertainties must 
be reduced about three orders of magnitude. Since 
the a priori uncertainties are of the order of one 
percent of the nominal thrust, this means that the 
instrumentation must be capable of measuring the 
low thrust accelerations (which are of the order of 
10“5 g) to an accuracy of 0. 001 percent. The 
alignment of the accelerometer axes in inertial 
space must be known to a comparable accuracy, or 
about 10“5 radian. These requirements appear to 
be well beyond the state of the art. 6) Further- 
more, there is the question of reliability. In each 
of two flight tests of a low-g electrostatic acceler- 
ometer, the instrumentation failed after a brief 
period of operation. (5,6) 

Fortunately, there is another approach that is 
much more promising, namely, the development 
of an Earth -based data type that is relatively 
insensitive to stochastic acceleration effects. 
Ondrasik and Rourke (7,8,9) have shown that dif- 
ferencing simultaneous range -rate data from widely 
separated stations produces a data type that con- 
tains the essential right ascension and declination 
information and is largely unaffected by stochastic 
accelerations. This technique is called "two- 
station tracking" or QVLBI (Quasi -Very Long Base- 
line Interferometry). Although still in its infancy 
and not yet proven in operations, this simple con- 
cept promises to be the breakthrough needed for 
low thrust navigation to become competitive with 
ballistic. 

With the emergence of QVLBI as a viable track- 
ing strategy a re-evaluation of low thrust orbit 
determination capability is needed, taking into 
account recent developments in SEP thrust sub- 
system error modeling. There are several 
important questions to be answered: 

1) How effective is QVLBI compared with con- 
ventional tracking for SEP missions? 

2) Is explicit differencing of the simultaneous 
data necessary to achieve the desired effect? 


Most of the material in this memorandum also appeared as AIAA Paper No. 73-174, "Earth- Based 
Orbit Determination for Solar Electric Spacecraft with Application to a Comet Encke Rendezvous," 
presented at the AIAA 11th Aerospace Sciences Meeting in Washington, D. C. , January 10-12, 1973 
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3) What degradation of performance occurs as 
the result of an incorrect thrust subsystem error 
model ? 

4) What are guidelines for conservative error 
modeling; i, e, , are some types of model errors 
more damaging than others? 


deviation of the bias, standard deviation of the 
time -varying component, and correlation time, 
respectively, for the thrust magnitude error; 

^<j)» '^(j) corresponding parameters for 

the thrust orientation angles. The two angles are 
assumed to have identical statistical properties so 
that a single set of parameters applies to both. 


5) What is the relative importance of the 
various error sources? 

This paper provides preliminary answers to 
these and related questions, using a 1980 rendez- 
vous with the Comet Encke as a representative 
SEP mission. First, the basic thrust subsystem 
stochastic error model is described and baseline 
values for the error model parameters are given. 
Five different tracking strategies, including 
QVL.BI, are defined and applied to a 15 day segment 
of the 950 day Encke mission, using the baseline 
values of thrust subsystem error model param- 
eters to determine the baseline orbit determination 
performance. Then, for three of the five tracking 
strategies, the error model parameters are varied. 
Both optimal and suboptimal performance are 
determined for a wide range of parameter values. 
The total orbit determination errors for QVLBI and 
single -station tracking are also broken down and 
analyzed according to individual error sources, 

II. Thrust Subsystem Error Model 

The results of any SEP orbit determination 
study depend very strongly upon the assumptions 
made about the stochastic nongravitational acceler- 
ations, For this reason, considerable effort has 
been directed by the SEP Navigation Development 
Team at JPL toward identifying the major thrust 
subsystem error sources and developing a satis- 
factory model for their representation. The result- 
ing “process noise model” is the basis for the 
present study. The basic characteristics of the 
error model will be outlined in this section. 


To obtain baseline values for the standard 
deviations of the process noise model it was neces- 
sary to identify the contributing error sources, 
determine the expected error in each parameter, 
and compute the corresponding effect on the thrust 
process. For the orientation angles the error 
analysis is relatively simple. The bias in the 
angles is due mainly to grid warpage and is typically 
about 2 degrees (, 035 radian). The time -varying 
angle components represent the random pointing 
error of the spacecraft within the deadband of the 
attitude control sensors. This is expected to be 
less than one degree, or approximately . 01 radian. 

Analysis of thrust magnitude errors is somewhat 
more complicated, as there are a number of thrust 
subsystem parameters contributing to the overall 
magnitude error. The nominal thrust magnitude is 
given by 

/ t ). + ^/T 7 }^\ 

T = ^r,|+2r,^ ) ^ ® ^ 

where 

It 3 = beam current 

Id 

Vg = net accelerating potential 
cos 9 = beam divergence factor 

= singly ionized fraction of mass flow 
^2 - doubly ionized fraction of mass flow 


The stochastic accelerations are assximed to be 
due solely to variations in the thrust vector, which 
is represented by its magnitude and two orientation 
angles (clock and cone angles with respect to the 
sun-probe vector). The errors in each orientation 
angle and in thrust magnitude are assvimed to be 
independent and to consist of a zero -mean time- 
varying random component superimposed upon a 
bias. Each time -varying component is modeled as 
a first order Gauss -Markov random process; (10) 
i, e. , each component is assumed to satisfy a 
stochastic differential equation of the form 


f = charge exchange parameter 

K = constant 

A linear expansion of Eq, (3) about nominal pa ram 
eter values gives the following expression for the 
relative thrust magnitude error (AT/T). 

th +^1 


AT 

T 


AI„ , AV„ 


^ - (1/t)^ + w (1) 

where w is a Gaussian white noise process such 
that 


E[w(t)] = 0 

Ejw(t^)w(t2)j = (Zo-^/t) S 


( 2 ) 


Each random process of this form is completely 
characterized by two parameters: the standard 
deviation (a), which is a measure of uncertainty, 
and the correlation time (t), which is a measure of 
transience. The following system of notation will 
be used; ot, and Tq» denote the standard 


+ ^ (4) 

cos 0 

From Eq. (4) and from available information 
about the range of parameter variations, the con- 
tribution of each parameter to the overall thrust 
magnitude error (as a percent of nominal thrust) 
can be determined. This information is sum- 
marized in Table 1. The dominant error sources 
are statistically independent, so the total thrust 
error is the root -sum -square (RSS) of the indivi- 
dual errors. 

The standard deviations in Table 1 are for a 
single thruster. When more than one thruster is 
operating, the standard deviations for the process 
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Table 1 Thrust Subsystem Error Summary 


Parameter 

Nominal 

Value 

Calibration 

Accuracy 

(%) 

Variation 
with Time 
(%) 

Contribution to Thrust Error 
(% of Nominal Thrust) 

"■Tb 

<Tt 

’b 

Programmed 

± 

0.5 

± 1. 5 

0. 5 

1. 5 


1500 V 

± 

0.5 

± 1.0 

0.25 

0.5 

cos 6 

0. 96 

± 

2.0 

± 3.0 

2.0 

3.0 


0.80 - 0, 85 

± 

1.0 

±5.0 

0. 02 - 0. 05 

0. 10 - 0.25 


0. 04 - 0. 07 

± 20.0 

± 25.0 

0.5 - 1.0 

0.50 - 1.25 

C 

0.995 

± 

0. 15 

± 0. 15 

0. 15 

0. 15 

Total (RSS) 




2.2 

3.5' 


noise model are as given in Table 2. The total 
error is the RSS of individual thruster errors in 
each case with the exception of the pointing error 
from the attitude control system (cJ^), which is 
independent of the number of thrusters. 

Table 2 Process Noise Model Standard Deviations 


Number of 
Operating 
Thrusters 

Thrust Magnitude 
(% of Nominal Thrust) 

Orientation Angles 
(radians) 

"■t 

*^Tb 


®^<|>b 

1 

3. 5 

2.2 

0. 010 

0.035 

2 

2. 5 

1. 6 

0. 010 

0. 025 

3 

2.0 

1. 3 

0. 010 

0. 020 

4 

1. 75 

1. 1 

0.010 

0.0175 

5 

1. 6 

1. 0 

0. 010 

0. 016 


The present study is based on a 15 day segment 
of the 950 day Encke rendezvous mission terminat- 
ing 5 days before encounter. During this time 
power matching requirements dictate that there 
will be four thrusters operating. Thus, the stan- 
dard deviations from Table 2 corresponding to four 
thrusters were used as baseline values. The indi- 
cated thrust magnitude error of 1. 75 percent 
represents acceleration errors of . 6 to , 8 x 10“® 
km/sec^ during this particular mission segment. 

The correlation times required for the error 
model are somewhat more nebulous than the stan- 
dard deviations. Observed variations in the quan- 
tities contributing to the thrust magnitude error 
indicate relatively long correlation times, of the 
order of days to weeks. Since the longer correla- 
tion times in this range approach a bias effect, of 
which the orbit determination process is relatively 
tolerant, in the interest of conservative modeling a 
value for T-p toward the lower end of this spectrum, 
namely, 5 days, was selected as the baseline value. 
Pointing errors are higher in frequency and more 
transient than the thrust magnitude fluctuations. 
Indeed, correlation times of the order of hours 
appear to be appropriate. Accordingly, a baseline 
value of 3 hours was selected for t^. 

Conspicuously absent from this error model is 
any mention of mass variations, which, if present, 
of course, also contribute to the stochastic acceler- 
ation. The fundamental thrust subsystem param- 
eters do not affect mass directly, but three of them, 
namely, I B# ^1 sind ultimately affect mass 
through variations in the mass flow rate m. Indeed, 
a linear expansion for m analogous to Eq. (4) indi- 
cates that the time -varying error component in m 
may be as large as ±5%. This seemingly large 


error is actually negligible for the following 
reasons: First, the variation in fuel (i.e., mass) 
expended over the time period of a typical tracking 
interval is only a small portion of the total fuel for 
the 950 day mission, which in turn is a relatively 
small portion (about 30%) of the total spacecraft 
mass; e. g. , for the particular 15 day period of 
interest in this study, a 5% variation in mass flow 
rate represents a maximum of . 09% variation in 
the total mass, which is clearly negligible with 
respect to the thrust magnitude errors. Further- 
more, since the mass variation is the integral of 
the error in m, higher frequency components are 
attenuated, and the slowly varying components, 
even if they were not of negligible magnitude, could 
be adeqxiately represented as a bias over a short 
tracking interval. The foregoing observations do 
not preclude the possibility of a large initial 
uncertainty in mass due to the cumulative effect of 
mass flow rate errors over a long period of time 
preceding the data arc. In this study an initial 
uncertainty of 100 kg (out of a total spacecraft 
mass of 1200 kg) was assumed. 

III. Filter Configuration and Tracking Strategies 

The orbit determination filter used for this 
study was a discrete form of the Kalman sequential 
filter, (2, 10) with the additional capability of ’’con- 
sidering'’ the effect of specified parameters with- 
out explicitly estimating them. A companion 
evaluation algorithm allowed separate specification 
of the process noise parameters for the filter model 
and the assumed real-world environment. The 
stochastic thrust errors were treated as piecewise 
constant functions over one hour time intervals, 
closely approximating the continuous Gauss -Markov 
process of Eq. (I ) in the manner suggested by 
Russell and Curkendall. (H) The details of the 
filter algorithm are given in the Appendix. 

The estimated parameters, in addition to the 
spacecraft state, were the initial mass, the biases 
and stochastic components of the two thrust vector 
orientation angles, and the stochastic component of 
the relative thrust magnitude AT/T. The bias in 
relative thrust magnitude was not included because, 
under the assvimption that the nominal mass is 
constant during the tracking interval (it actually 
decreases about 2% in this case), the effect of a 
bias in relative thrust magnitude is indistinguishable 
from that of a bias in the mass (i. e. , they are per- 
fectly correlated). Therefore, it is redundant for 
the filter to explicitly estimate both quantities. 

The large initial mass uncertainty (100 kg) assumed 
for this study effectively absorbs the bias compo- 
nent in thrust magnitude, allowing the latter to be 
neglected. A priori uncertainties for all the above- 
mentioned quantities are given in Table 3, along 
with the assumed data accuracies. The DSN sta- 
tion locations were considered as an error source, 
but not explicitly estimated by the filter. 

Two single -station (SS) and three multi -station 
(MS) tracking strategies were investigated in this 
study. For ease of reference in the sequel, the 
following abbreviated terminology will be used to 
refer to specific combinations of data types. 

1) SS Doppler -only: two-way range -rate data 
from a single station. 

2) SS conventional: two-way range and range - 
rate from a single station. 
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Table 3 Key Assumptions for the Baseline Case 


Description 

Numerical Value 

Initial S/C State Uncertainty 


RSS Position 

17, 300 km 

RSS Velocity 

17. 3 m/s 

Mass 

100 kg 

Thrust Magnitude Errors 


Bias Component 

(absorbed in mass) 

Time “Varying Component 

1. 75% 

Correlation Time 

5 days 

Orientation Angle Errors 


Bias Component 

0. 0175 rad 

Time- Varying Component 

0.010 rad 

Correlation Time 

3 hr 

DSN Station Location Errors 


Spin Radius 

1.5 m 

Longitude 

3. 0 m 

Parallel to Spin Axis 

5.0 m 

Data Noise 


Range Rate (for 1 min sample) 

1 mm/s 

Range 

3 m 

Differenced Range Rate 

\TZ mm/s 

Differenced Range 

m 


3) MS conventional: two-way range and range - 
rate from three stations; no simultaneous data, 

4) MS simultaneous: two-way range and range - 
rate from three stations plus simultaneous range 
and range -rate from pairs of stations during over- 
lapping view periods. 

5) MS differenced (or QVLBI): two-way range 
and range -rate from three stations plus explicitly 
differenced (two-way minus three-way) simul- 
taneous range and range -rate from pairs of stations 
during overlapping view periods. 

The analytical basis for the information content 
of MS differenced data is given in References 7-9. 

It should be emphasized that the latter two strate- 
gies use exactly the same data, but in essentially 
different ways. In the case of QVLBI the explicitly 
differenced data is weighted conventionally accord- 
ing to the data accuracies given in Table 3, but 
the two-way data in this case is de -weighted; i. e. , 
it is treated by the filter as if it were less accurate 
than it really is. The purpose of de-weighting 
is t o minimize the sensitivity to process noise 
inherent in the two-way data. The two-way data 
cannot be eliminated completely, as it provides 
the geocentric range and range -rate information 
that is largely lost from differencing the simul- 
taneous data. In this study the two-way range and 
range -rate data for QVLBI were weighted at 1 km 
and 250 mm/s (for a one minute sample), 
respectively. 

All single -station tracking was assumed to be 
done from Goldstone (DSS 14). Tracking stations 
at Madrid (DSS 61) and Johannesburg (DSS 51) were 
added for the multi -station strategies. These sta- 
tions were selected because of their favorably over- 
lapping view periods to provide the simultaneous 


two-way and three-way data required for QVLBI, 
as illustrated in Fig. 1. The basic tracking 
pattern was repeated once each day. 



Fig, 1 Simulated Multi -Station Tracking Pattern 
IV, Baseline Performance 

The baseline orbit determination results are 
given in Fig. 2. The RSS position error at the 
end of the data arc (E-5^) is shown for each of 
the five tracking strategies for three cases. The 
first case is the "optimal*' sequential filter (optimal 
in the sense that the stochastic thrust errors are 
correctly modeled by the filter*s process noise 
model). This set of bars represents the limiting 
orbit determination capability under the assump- 
tions of this study with the process noise as 
represented by the baseline model. There is 
about an order of magnitude improvement with 
each change of tracking strategy as one proceeds 
from SS Doppler -only to SS conventional, MS 
conventional, and MS simultaneous data, respec- 
tively. Since the process noise model is correct, 
implicit differencing of the simultaneous data by 
the filter yields good results, and nothing is 
gained by providing the filter with explicitly dif- 
ferenced data. Thus, the last two strategies give 
about the same performance in this case. 

To fully appreciate the power of QVLBI, one 
must look at something other than optimal filter 
performance. That is the main purpose of the 
second set of bars in Fig, 2. These represent 
the error in the state estimate resulting from a 
batch filter solution in the presence of SEP stochas- 
tic accelerations; i. e, , these are the errors that 
would result if the real-world stochastic accelera- 
tions were actually as represented by our baseline 
model, but the filter attempted to estimate all 
thrust variations (including the time-varying 
components) as if they were biases. (This is 
equivalent to using infinite correlation times in 
the process noise model. ) It is not intended to 
suggest that one would ever attempt orbit determina- 
tion for an SEP mission using a batch filter, but 
this is a simple means of demonstrating the rela- 
tive sensitivity of the various tracking strategies 
to incorrectly modeled (or unmodeled) stochastic 
accelerations. We have assvimed a very bad model, 
namely, a simple bias. In this case QVLBI (MS 
differenced data) is clearly superior, giving errors 
nearly two orders of magnitude smaller than any 
of the other tracking strategies (but still severely 
degraded from the optimal filter performance). 
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The MS simultaneous data in this case offers no 
advantage over conventional tracking. 


Process Noise Parameters 


V 



SEP WITH SEP WITH BALLISTIC 

SEQUENTIAL BATCH 

FILTER FILTER 


Fig. 2 Baseline Orbit Determination 
Pe r f o rmanc e 

As a basis for comparison, the ballistic orbit 
determination performance is also included in 
Fig. 2. These are the errors that would result if 
the Encke mission could be flown by a ballistic 
spacecraft following the SEP trajectory. We 
assume that the stochastic accelerations are negli- 
gible in this case so that data noise and station 
location uncertainties are the primary error 
sources. The orbit determination errors in this 
case are qujte small even by ballistic standards 
due to the favorable geometry of the trajectory 
near Encke encounter. (The spacecraft is about 
0. 3 A.U, from Earth with a declination of 
60 degrees. ) Also, there is little basis for choos- 
ing one tracking strategy over another in this case. 
Even SS Doppler -only tracking gives good perfor- 
mance. The presence of process noise severely 
degrades the filter *s ability to infer position infor- 
mation from range -rate data. This accounts for 
the dramatic disparity between the ballistic and 
SEP cases for Doppler -only tracking. Only in the 
case of optimal sequential filtering with MS simul- 
taneous or MS differenced data are orbit determina- 
tion errors for SEP within an order of magnitude 
of those for ballistic spacecraft. 


In this section we consider the effect of changes 
in the process noise standard deviations and corre- 
lation times, but always under the assumption that 
the filter *s model is a correct representation of the 
actual stochastic environment, thus the optimal 
orbit determination capability is determined as a 
function of process noise parameters. For this 
part of the investigation three tracking strategies 
were considered; SS conventional, MS simultaneous, 
and MS differenced. The correlation times were 
held fixed at baseline values as the standard devia- 
tions were varied, and conversely. 

First, consider the sensitivity to standard devia- 
tions as given in Fig. 3, The standard deviations 
here are expressed as a ratio with respect to base- 
line values; i. e. , the three process noise standard 
deviations (two angles and relative thrust magni- 
tude) were always varied by the same multiplica- 
tive factor. As expected, the simultaneous and 
differenced data give essentially the same per- 
formance. Of course, the optimal performance of 
differenced data is always slightly inferior to that 
of simultaneous data because of the de -weighting of 
two-way data in connection with the former. Note, 
however, that as the process noise level increases, 
the separation of the two curves decreases. 

The optimal performance for SS conventional 
tracking is not only orders of magnitude inferior to 
the MS strategies, but is also much more sensitive 
to the process noise level; i. e, , the slope of the 
curve is much greater for SS conventional data. 

Since the curve in this case is nearly a straight 
line on the log plot, it can be fit fairly well with an 
empirical formula of the form 

= k(<r/(T^P (5) 

where 

= RSS position error, km 

' ratio of process noise standard 
deviations to baseline values 

k = 3000 for SS conventional data 

^ = 0. 89 for SS conventional data 

The value of the exponent r\ being near unity indi- 
cates that orbit determination errors are almost 
linear with the process noise level. The empirical 
formula (5) does not fit the QVLBI case as well, 
but would be applicable in the neighborhood of the 
baseline values with k = 55 and rj = 0, 22. 

The sensitivity to correlation times is given in 
Fig, 4. For this part of the investigation the three 
process noise correlation times were assumed to 
be the same, but a sufficiently wide range of varia- 
tion was considered to include the baseline values 
of both the orientation angles (0. 125 day) and thrust 
magnitude (5 days). Once again, SS conventional 
tracking exhibits the greatest sensitivity. The per- 
formance of MS simultaneous and MS differenced 
data is relatively flat with respect to correlation 
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RSS POSITION ERROR, km 



Fig. 3 Sensitivity of Optimal Orbit Deter- 
mination Performance to Process 
Noise Standard Deviation 


time with the maximum error occurring near one 
day. Conventional tracking, on the other hand, 
gives a relatively sharp peak at 0. 2 day (about 
5 hours) and rapid improvement for larger cor- 
relation times (approaching a bias). 


VI. Sensitivity to Mismodeled 
Standard Deviation 


In this section and the next we consider the 
more realistic situation where the actual stochas- 
tic environment differs from the assumed process 
noise model of the filter. This kind of analysis 
not only gives a more realistic assessment of 
orbit determination capability, but also leads to 
guidelines for conservative modeling. 

Figures 5 and 6 give the orbit determination 
errors for conventional data and QVLBI (MS 
differenced data), respectively, for a wide range 
of actual and assumed standard deviations. The 



CORRELATION TIME, days 

Fig. 4 Sensitivity of Optimal Orbit Deter- 
mination Performance to Process 
Noise Correlation Time 

dashed line in each of these figures is the locus of 
minima of the family of curves and is identical with 
the corresponding optimal performance curve from 
Fig, 3. In the case of conventional data (Fig. 5) 
the dashed optimal performance curve divides the 
sub optimal performance into two regimes. To the 
left the curves are relatively flat, indicating that 
using a standard deviation that is too low in the fil- 
ter's model does not degrade performance signifi- 
cantly from the optimal. Using too large a stan- 
dard deviation, on the other hand, can degrade per- 
formance severely, as indicated by the rapid rise of 
curves to the right. This pattern does not hold for 
QVLBI, however. In this case (Fig. 6), the sub- 
optimal performance curves tend to be more nearly 
symmetric about their minima. 

The suboptimal performance of MS simultaneous 
data with respect to mismodeled standard deviation 
was found to be virtually indistinguishable from that 
of MS differenced data as. given in Fig. 6, and so 
those results are not repeated here. At one 
extreme, with org^/o'b “ ^ind = 0. 1, the 

orbit determination errors were 6 percent larger 
with simultaneous data. At the other extreme, with 
cTa/o-b = 0. 25 and = 10, the errors were 

5 percent smaller. Thus, if the assumed standard 
deviation is more than an order of magnitude too 
low, explicit differencing may offer a very slight 
advantage. 

The suboptimal orbit determination performance 
can be better understood by analyzing the individual 
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<7^ = ACTUAL STD DEV 
(Tf = FILTER ASSUMED STD DEV 
= BASELINE VALUE OF STD DEV. 


/ a-h = 4 




TOTAL ( <T„/ 0-. = 1) 



Sensitivity of SS Conventional Data to 
Mismodeled Standard Deviation 


Fig. 7 Contribution of Individual Error 
Sources to Total Orbit Determination 
Error for SS Conventional Data 


O-^ = ACTUAL STD DEV 

o-f = FILTER ASSUMED STD DEV 

0-. - BASELINE VALUE OF STD DEV 




noise. Station locations and data noise are 
negligible error sources in comparison. As the 
assumed standard deviation is increased, the fil- 
ter does a better job of estimating state deviations 
due to the stochastic thrust components, but the 
error due to the bias components increases rap- 
idly, with consequent overall degradation. Note 
that the a priori state uncertainties are still sig- 
nificant even after 15 days of tracking. Indeed, if 
the assTimed standard deviation is much too large, 
the a priori uncertainties are dominant, indi- 
cating that the filter is getting very little informa- 
tion from the data. 


Fig. 6 Sensitivity of MS Differenced Data to 
Mismodeled Standard Deviation 

error sources. This is done in Fig, 7 for con- 
ventional data and in Fig. 8 for QVLrBI. Fig. 7 
shows that with conventional data orbit determina- 
tion errors are dominated by the effect of process 


In contrast, from Fig, 8 QVLBI is found to be 
data noise limited, in general. Process noise 
becomes the dominant error source only when the 
assumed standard deviation is much too small. 
Station locations are a significant, but never domi- 
nant, error source, and their effect is relatively 
insensitive to the process noise model. The effect 
of a priori state uncertainties after 15 days of 
tracking is too small to appear in Fig. 8. 

VII. Sensitivity to Mismodeled Correlation Time 

The preceding section dealt with the case 
where the dynamic structure of the process noise 
model (as reflected by the correlation time) was 
correct, but the level of process noise (as 
reflected by the standard deviation) was possibly 
incorrect. In this section we consider the opposite 
circumstance. This is perhaps the most realistic 
form of mismodeling, as we usually know more 
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Fig. 8 Contribution of Individual Error 
Sources to Total Orbit Determination 
Error for MS Differenced 
Data (QVLBI) 

about the magnitude of stochastic forces than 
about their underlying dynamics. The results are 
displayed in Figs. 9-11. Orbit determination 
errors here are plotted against actual correlation 
time for various assumed values (cf. Figs, 5 and 6). 
The primary objective, of course, is to see if any 
assTimed value gives satisfactory performance over 
a wide range of actual values. As before, to sim- 
plify the analysis in this part of the investigation, 
we have disregarded the difference between orien- 
tation angle and thrust magnitude correlation times 
in the baseline thrust model and have assumed all 
three components of the process noise to have the 
same correlation time, both in the filter model and 
in the supposed actual environment. 

The results for conventional data are given in 
Fig, 9. If the lower envelope for this family of 
curves were constructed, it would be precisely the 
optimal performance as a fiinction of correlation 
time given earlier in Fig. 4. The spread of curves 
above the lower lenvelope indicates the extent of 
potential deterioration of performance from the opti- 
mal, assuming an incorrect correlation time. The 
curves are very widely spaced to the right, indicat- 
ing considerable loss of performance if a short cor- 
relation time is assumed when the actual is large. 
The spread is much smaller to the left; however, 
the optimal performance is much worse here; 
therefore, less deterioration can be tolerated. If 
the actual correlation time is completely unknown, 
a conservative choice for the filter model appears 
to be 0.2 day, as this choice makes the orbit 



0.01 0.1 1 10 100 
days 


Fig, 9. Sensitivity of SS Conventional Data to 
Mismodeled Correlation Time 



Fig, 10. Sensitivity of MS Differenced Data to 
Mismodeled Correlation Time 
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Fig. 11. Comparison of MS Simultaneous and MS Differenced Data with 
Mismodeled Correlation Time 


determination process least sensitive to the actual 
correlation time. Furthermore, the peak error 
with Tf = 0. 2 occurs at Ta = 0. 2, This means 
that if pre -mission analysis is done by computing 
optimal performance with rf = 0, 2, the actual 
performance can be no worse regardless of the 
actual correlation time. 

The analogous family of curves for QVLBI is 
given in Fig. 10, where the lower envelope indicat- 
ing optimal performance is shown as a dashed 
line. A correlation time of one day appears to be 
an excellent choice for the filter model in this 
case. The curve corresponding to rf = 1 day is 
remarkably flat and very close to the optimal per- 
formance curve over the entire range of actual 


correlation times. One question remains: How 
does simultaneous data compare with explicitly 
differenced data when the correlation time is mis- 
modeled? The answer is given in Fig. 11 where 
the performance of these two strategies is com- 
pared for four different values of the assumed cor- 
relation time. The explicitly differenced data is 
found to be superior only when the actual correla- 
tion time is at least an order of magnitude smaller 
than the filter *s assumed value. For a very large 
assumed correlation time (e. g. , 100 days), the 
improvement with explicitly differenced data may 
be considerable (as much as an order of magnitude), 
but the smaller the assumed correlation time, the 
less difference there is between the two strategies. 
For Tf = 1 day they give comparable performance. 
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VIII. Concluding Remarks 


References (contd) 


New multi- station tracking techniques have 
been shown to be superior to conventional single- 
station tracking for SEP missions, reducing orbit 
determination errors by orders of magnitude and 
making SEP navigation capability competitive with 
ballistic. Explicitly differenced multi- station 
data is least sensitive to modeling errors and may 
give satisfactory performance even when the 
process noise is grossly mismodeled. However, 
if a reasonably good process noise model is avail- 
able, or if a conservative model is chosen (i. e. , 
a relatively short correlation time is assumed), 
explicit differencing of the two-way and three-way 
data is not required. Although much larger 
errors result from single- station conventional 
data, this tracking mode undoubtedly will be used 
during the cruise portion of an extended SEP 
mission, where errors are less critical, to mini- 
mize the cost of operating the DSN. Therefore, 
the guidelines developed herein for effective use 
of conventional data may also be of more than 
academic interest. 
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APPENDIX 


The Filter Algorithm 


The filter algorithm assumes a linear dynam- 
ic system of the form 



k-1 


(A-1) 


with scalar observations 


Zk = [ A 



(A-2) 


where is the vector to be estimated, yj^ is the 
vector of uncertain parameters and stochastic 
variables to be "considered," but not explicitly 
estimated, and wj^ and Vj^ are white noise proc- 
esses such that 


E[wj^] = = 0 


(A-3) 


Equations (A-5) - (A-8) are intended to show the 
essential characteristics of the algorithm as 
simply as possible. Of course, more efficient 
means of implementation are used in practice. 

This algorithm is essentially the Schmidt- Kalman 
filter 2, 1 3 ) with two minor differences. First, 
the dynamic model of Eq. (A-1) is slightly more 
general, allowing the considered variables to be 
dynamic or even stochastic rather than simple 
biases. Second, the term A Pj^A^ in Eq. (A-6) is 
replaced by H P'hT in the Schmidt-Kalman filter. 
This difference is probably not significant in prac- 
tice, but Eq. (A-6) represents the more conser- 
vative choice, as the denominator of the filter 
gain is smaller in this case, giving a larger value 
of K and consequently a little more protection 
against filter divergence. 


Now suppose the covariances P(0), Q, and R 
input to the filter do not represent the true sta- 
tistics of the ensemble of random variables in- 
volved. but r^her the true statistics are given by 
P(0), Q, and R, respective^. Then the propaga- 
tion of the true covariance P(k) will be given by 


where the dependence of U, V, M, A, B, Q, and 
R upon the index variable k has been suppressed 
to simplify the notation. Let Xj^ be the estimate 
of xj^ and define the error covariance 


P(k) = E 



(k) 

pT 

L xy(k) 


Py(k) 


(A-4) 


where 


X - X 



Given xq and P(0), the recursions for propagating 
the estimate and covariance are 


P'(k) = $P(k - 1)$'^ + Q (A-5) 

K= [P]^(k) PJ^^(k)] h'^[A P^(k)A'^ + R (A-6) 

P(k) = [I - GH]P'(k)[l - GH]^ + GRG^ (A-7) 


where 


$ = 


V 

M 


H = [A B] 


G 


K 

0 


(A-9) 


P'(k) = $P(k - 1)$'^ + Q (A-10) 

P(k) = [I - GH]P'(k)[l - GH]'^ + GRg'^ (A-11) 


where the filter gain K (and therefore G) is as 
computed from Eq, (A-6). Thus, the addition of 
Eqs. (A-10) and (A-11 ) to the algorithm allows 
parallel propagation of covariances for the "filter 
world" and the "real world," as required for this 
study. 


It may not be apparent that the simple scheme 
described above allows for evaluation of mis- 
modeled dynamics as well as mismodeled sta- 
tistics. Suppose, for example, that the effect of 
mismodeled correlation time for the random com- 
ponent of thrust magnitude is to be evaluated. 

Then the thrust magnitude dynamics would be in- 
cluded redundantly in the linear system model 
(i. e. , the the V and M matrices) using both cor- 
relation times. The a priori and driving white 
noise variances corresponding to one correlation 
time would be set to zero in P(0) and Q, and those 
corresponding to the otherjc or relation time would 
be set to zero in P(0) and Q. As a result, the 
filter gain is based on one dynamic model, and 
the propagation of the true covariance takes place 
with a different dynamic model. Thus, the dynam- 
ic modeling error is appropriately accounted for. 
This is the technique used for the present study. 
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